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1.  INTRODUCTION 


In  a  recent  study,  ^  a  formalism  and  computer  code  (EMA.BIC)  were 
developed  to  invert  profiles  of  plume  infrared  radiance  and  absorptance  obtained 
in  transverse  scans  across  the  lateral  extent  of  test  engine  plumes  near  the 
nozzle  exit  plane  in  order  to  retrieve  the  radial  profiles  of  temperature  and  gas 
species  concentration  in  the  plume  flow  field.  Principal  aspects  of  the  formula¬ 
tion  include:  (1)  the  use  of  state-of-the-art  band  radiation  models  that  treat 
highly  nonuniform  plumes  and  account  for  Doppler  broadening  effects;  (2)  the 
development  of  an  iterative  Abel  inversion  algorithm  that  allows  the  application 
of  the  well-known  Abel  inversion  to  nonthin  optical  sources;  and  (3)  the  formula¬ 
tion  of  a  data  smoothing  procedure. 

Consideration  of  smoothing  is  necessary  to  a  successful  inversion 
algorithm.  Regardless  of  the  inversion  method  used,  experimental  random 
noise  on  the  input  emission/absorption  (E/A)  profiles  is  amplified  in  the  inversion 
process  and  can  produce  radial  profiles  so  noisy  that  they  are  useless.  The 
smoothing  procedure  used  in  EMABIC  is  a  straightforward  presmoothing  of  the 
input  transverse  E/A  data  before  inversion. 

The  occurrence  of  random  experimental  error  on  the  input  profiles  implies 
an  associated  error  on  the  retrieved  radial  profiles.  In  the  previous  work,  the 
propagation  of  errors  through  the  inversion  process  was  studied  empirically 
by  superimposing  random  errors  on  smooth  transverse  profiles  and  analyzing 
their  effect  on  the  retrieved  radial  profiles.  In  the  present  work,  the  formal 
theory  of  linear  error  propagation  by  transformation  of  the  variance-covariance 


1.  S.  J.  Young,  Inversion  of  Plume  Radiance  and  Absorptance  Data  for 

Temperature  and  Concentration,  AFRPL-TR-78-&0,  Air  Force  Rocket 
Propulsion  Laboratory,  Edwards  Air  Force  Base,  Calif.  29  September  1978. 
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error  matrix  of  the  data  is  applied.  The  resultant  error  propagation  algorithm 
has  been  coded  into  EMABIC  so  that  automatic  error  analysis  can  be  made  along 
with  inversion. 

(2) 

The  general  theory  of  linear  error  propagation  is  well  established;  ’  and 

only  an  outline  and  final  results  are  presented  here.  This  outline  and 

application  of' the  results  to  the  present  problem  are  treated  in  Section  2.  A 

(3  4) 

similar  treatment  of  this  problem  has  been  made  by  Limbaugh.  v  '  '  An 
example  application  of  the  method  and  comparison  of  results  with  the  previous 
empirical  results  are  made  in  Section  3.  Preparation  of  control  cards  to  operate 
EMABIC  in  the  error  propagation  mode  is  described  in  an  appendix. 


2.  A.  A.  Clifford,  Multivariate  Error  Analysis,  Applied  Science  Publishers 
Ltd.,  London,  1973. 

3.  C.  C.  Limbaugh,  An  Uncertainty  Propagation  Analysis  for  an  Infrared 
Band  Model  Technique  for  Combustion  Gas  Diagnostics,  AEDC-TR-76- 155, 
Arnold  Engineering  Development  Center,  Arnold  Air  Force  Station,  Tenn., 
April  1977. 

4.  R.  T.  Shelby  and  C.  C.  Limbaugh,  Smoothing  Technique  and  Variance 
Propagation  for  Abel  Inversion  of  Scattered  Data,  AEDC-TR-76  - 163, 
Arnold  Engineering  Development  Center,  Arnold  Air  Force  Station,  Tenn., 
April  1977. 
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2.  ERROR  PROPAGATION  THEORY 


2.  1  General  Theory 

Consider  a  set  of  n  equations  defining  n  parameters  in  terms  of 
n  observables  X- 


yi  *x: 


Xn> 


=  1, 


(1) 


The  Xj  represent  measurable  (observable)  quantities, and  we  assume  that  we 
have  a  complete  description  of  the  measurements  in  the  form  of  a  set  of  mean 
values  x.(i=l,n),  standard  deviations  a(x-){i=l,n)  and  correlation  coefficients 
p(Xj*  Xj)  (i,  j  =  l.n).  This  last  quantity  is  a  measure  of  the  degree  of  dependence 
between  x.  and  Xj.  For  i  =  j,  p(xr  Xj)  s  1. 

The  mean  values  of  the  parameters  are  calculated  from  Eq.  (1)  with  the 
Xj  used  in  place  of  the  Xj-  Thus, 

y-  =  Yi(X1,  X2,...,Xn)  i  =  l.n  (2) 

The  functional  dependence  of  the  y^  on  the  x^  can  be  quite  general  for 

computing  the  mean  values  y..  In  order  to  treat  error  propagation,  however, 
we  invoke  a  linearity  approximation  and  write  Eq.  (1)  to  first  order  as 


n 


8yi  =£ 


iyi 


J-1  ^  ^ 


i  =  1 ,  n 


(3) 


where  the  partial  derivatives  are  evaluated  at  the  x .  ,  6y.  =  y.- y^  and  &Xj  =Xj  -  X  j. 
For  this  approximation  to  be  reasonable,  either  the  6 x j  rnust  be  small  or  the 
equations  must  be  nearly  linear  (so  that  second  derivatives  are  small).  Now 
define  the  partial  derivative  matrix  A 


r 


A  = 


1 


5*1 


ax. 


ay, 

ax. 


&y2 

ax. 


^1 


ay^ 

ax. 


(4) 


and  the  error  matrix  of  the  observables 


M  = 
X 


a2(Xj) 


o(Xj)  ct(X2)  p(X,»  X2)  •  •  •  ^(Xj)  d(Xn)  p(X  | ,  Xn) 


ct(X2)  cr(Xj)  p(x2»  xx)  o  (X2) 


CT(Xn)  CJ"(x 2 )  P(X^»  Xj) 


(5) 


This  matrix  is  the  variance-convariance  matrix  for  the  observable  array 
Xj  [ct  (X^)  is  the  variance  of  X.  and  o(X.)  cr(Xj)  p(X^,  Xj)  is  the  covariance 
between  X.and  Xj  ]  .  The  matrix  is  symmetric  since  p(Xj,  Xj)  =  p(Xj,  X.). 

The  corresponding  variance -covariance  matrix  for  the  parameter  array 

Yi  is 


M  =  AM  AT 

y  x 


(6) 
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V 


—  t  —  — 

where  A  is  the  transpose  of  A.  Explicitly,  the  ij  element  of  is 


n 

=  E  ^ 

k=l 


CJ(XX)  Eli  p(xk’  xi> 
1  =  1  *xz 


(7) 


with  Xk)  =  1  •  The  matrix  M  has  exactly  the  same  structure  as 

M  .  The  standard  deviation  of  y.  is 
X  yi 


°<yi)  =  V*’Vii 


(8) 


and  the  correlation  coefficient  connecting  y^  and  y.  is 


(M  ).. 

0(V  V  )  -  _ OI_ 

P  Yi’  Vj  cr(y.)  a(Yj) 


(9) 


In  the  following  sections,  this  general  result  will  be  applied  to  the 
propagation  of  errors  through  the  data  presmoothing  operation  and  then 
through  the  profile  inversion  operation  used  in  EMABIC.  In  the  first 
application,  the  variance -covariance  matrix  is  that  of  the  raw, 
unsmoothed  E/A  data,  and  is  the  corresponding  matrix  for  the  smoothed 
E/A  data.  In  the  second  application,  of  the  smoothed  data  becomes  the 

input  matrix  M  ,  and  the  computed  M  describes  the  error  characteristics 

x  v 

of  the  recovered  temperature  and  concentration  profiles  of  the  plume. 
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2.2  Error  Propagation  through  Presmoothing 


Consider  a  function  g(z)  defined  on  the  grid  z^(i  =  l,  M; 

0  =  Zj  <  z2  <  .  .  .  <  z^  =  R;  R  =  plume  radius)  by  g.  =  g(z,.).  The  array 
g.  represents  either  the  unsmoothed  transverse  data  sets  IsL  (radiance) 
or  1  -  T\  (absorptance).  Let  f.  be  the  smoothed  version  of  g.  (and  defined 
on  the  same  grid).  The  smoothing  routine  of  EMABIC  defines  the  f. 

f  =  Hg'  (10) 


where  g1  is  the  array  gj  =  d.  g.  (i  =  1 ,  M  -  1),  and  H  is  the  inverse  of  the 
(M  -  1)  X  (M  -  1)  banded  symmetric  matrix 


kl 

k2 

k3 

k4 

k5 

k6 

k7 

7? 

00 

k9 

k10 

kll 

k12 

kl  3 

k14 

k15 

kl6 

k  = 

00 

•  •  • 

k5M-28 


(ID 


•  •  • 

k5M-23 

k5M-22  k5M-21 

k5M-20 

k5M-l  9 

k5M-18 

k5M-l  7 

k5M-16 

k5M-15 

k5M-14 

k5M- 1  3 

k5M-12 

k5M-ll 

i 

J 

! 


The  elements  of  K  are 


kl  =  Z+  bl^  Y  +  dl 
k2  =-<a2b2  +  2bl  C1)Y 


k4  "  k2 

k5  =  (a  32+  b22+  4Cj2)Y  +  d2 
k6  =“(a3b3  +  b2C2)  V 


*7 

3  , 

k5i-7  = 

k5i- 1 3 

k5i-6  = 

k5i-9 

k5i-5  = 

(a2+l  +b.2  +  ci21)7+d. 

k5i-4  = 

-<ai+l  b.  +  1  +b.c.)Y 

k_ .  ,  = 
5i-3 

ai+lCi+lY 

k5M- 1  3 

=  k5M- 1 8 

k5M- 1 2 

=  k5M  -14 
..2  2 

k5M- 1 1  =  (b  M-l  +  C  M-2*  y  +  dM- 
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where 


(with  zQ«-z2),  d1  =  Z2~ZV  di  =  <zi+1  -Zj.j)  (i=2,  M-l),  and  y  is  a  smoothing 
parameter.  (The  pass  through  the  loop  with  i  =  3  in  Eq.  (12)  will  incorrectly 
define  kg=k2.  This  must  be  corrected  to  kg  =  k^  after  completing  the  loop. 

Also,  the  last  pass  through  the  loop  will  define  ^5^,13  incorrectly,  but  this 
is  corrected  automatically  by  the  next  definition  after  the  loop. )  Note  that  the 
order  of  the  system  has  been  reduced  from  M  to  M-l  because  fM  is  constrain¬ 
ed  to  f.  ,  =  0  . 

M 

The  smoothing  procedure  is  derived  and  discussed  in  Ref.  1.  Briefly, 
the  method  consists  of  finding  the  function  f(z)  that  has  the  smoothest  overall 
curvature  but  which  does  not  differ  from  the  unsmoothed  function  g(z)  by  more 
than  the  estimated  error  of  g(z).  As  the  smoothing  parameter  varies  from 
y  =  0  to  00  ,  the  degree  of  smoothing  varies  from  none  to  infinite,  y  is  selected 
as  that  value  that  satisfies 


R 


R 

J  [f(z,  y)  -  g(z)]2  dr 
0 


2 

e 


where  e  is  a  measure  of  the  errors  in  g(z).  Smoothing  is  applied  separately  to 
N.  and  1-Tj  so  that  there  are,  in  general,  two  smoothing  parameters,  and 

yt  . 
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Equation  (10)  is  linear,  and  the  partial  derivative  matrix,  Eq.  (4), 

to  be  used  in  Eqs.  (6)  or  (7)  to  propagate  errors  to  the  smoothed 

profiles  is  just  the  matrix  H  =  K~*.  The  inversion  required  to  get  H  is 

(5) 

performed  in  EMABIC  with  a  standard'  '  Gaussian  elimination  routine 
designed  for  banded  matrices.  The  variance-covariance  matrix  for  the 
unsmoothed  profile  gJ  is  taken  as  the  diagonal  matrix 


J2  2  .  , 
djCT  (gx) 


d2*2(g2) 


,2  2/  \ 
d  M-1CT 


The  lack  of  nonzero  elements  off  the  main  diagonal  results  from  the  assumption 
that  the  transverse  data  have  been  collected  with  a  spatial  resolution  Az  that  is 
less  than  or  equal  to  the  data  step  size  across  OsSzfiR.  That  is,  Az  s  min 
(z^-Zj  j).  Then,  there  can  be  no  correlation  between  measured  values  of  g. 


and  g . , 


The  resultant  variance-covariance  matrices  for  the  smoothed  E/A  profiles 


=  H(yjj)  M^,H(y-)T 


M-  =  H(y?)  M-,  H(y?) ' 


System  1360  Scientific  Subroutine  Package  (360A-CM-Q3X)  Version  III, 
Programmers  Manual,  H20-0205-3,  IBM  Technical  Publication  Depart¬ 
ment,  White  Plains,  N.  Y.  ,  1968. 
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or,  explicitly, 


M-l 


(MN  >ij 
s  J 


=  X  hik(YN)G  (Nk}  hjk(YN)  dk 


k=  1 


M-l 


<MTs>ij  =  X  kik(v~)a2(V  hjk(v7>  d 


k  =  1 


(15) 


where  and  hjj(Y;p)  are  elements  of  H(y^)  and  H(y— ),  respectively.  Note 

that,  although  the  variance -covariance  matrices  of  the  unsmoothed  profiles  are 
diagonal,  the  matrices  for  the  smoothed  profiles  need  not  be.  With  increasing  y, 
the  correlation  between  elements  f.  and  f.  can  become  quite  large. 

In  order  to  handle  data  from  some  experimental  programs,  ^  EMAB1C 
allows  values  of  M  up  to  301.  Thus  H  could  be  as  large  as  300  x  300.  Both 
Me:  and  M—  could  presumably  then  also  be  each  300  x  300.  The  need  to 

lN  8  T  S 

retain  three  matrices  of  this  size  strains  the  storage  capacity  of  many  computers. 

Accordingly,  the  order  of  Me:  and  M—  is  reduced  in  size  to  NxN  where  N 

bl  s  Ts 

is  the  number  of  radial/transverse  zones  used  in  inversion  and  has  a  maximum 

value  of  50.  The  reduction  in  size  is  accomplished  by  simple  linear  interpolation 

on  the  four  elements  <Mfg).  (Mfg).  J+1,  (M£g)i+1>  .  and  (Mfg).  +  1>  j+1  that 

immediately  surround  the  grid  points  z^  =  (n-l)A  +  A/2,  z^  =  (m-l)  A  +  A/2 

and  M—  then 
T  s 

represent  the  variance -covariance  matrices  for  the  smootTied  profiles  at  the 
centers  of  the  N  transverse  zones  and  are  the  input  error  matrices  for  the 
subsequent  error  propagation  through  the  inversion  step. 


where  n,m=l,N  and  A  =  R/N.  The  reduced  matrices  M^j 
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2.  3  Error  Propagation  through  Inversion 


The  inversion  process  of  EMABIC  is  an  iterative  solution  of  the  equations 
of  radiative  transfer  defining  N(z)  and  t(z)  as  integrals  over  functions  of  T(r) 
and  c(r).  The  numerical  quadrature  approximations  to  these  equations  define 
N.  and  T.(i  =  1,  N)  in  terms  of  T^  and  cj(j-  1*  N)  by 


Ni  Ni  (Ti’  Ti+r***'  TN'  ci’  ci+l'  *  *  *  ’  CN) 
Ti  =  Ti  (Ti»  Ti+1  *  *  •  •  »  Tn»  ci»  ci+i»“*»  CN) 


i  =  1,  N  (16) 


The  inversion  process  solves  this  set  of  nonlinear  equations  to  give 


Tj  Tj  (Nj»  Nn,  t y  Tj+1>...,  tn) 

Cj  =  Cj  (*V  Tj'  Tj+1“--'  TN) 


j  =  1,  N  (17) 


From  Eq.  (3),  the  equations  describing  linear  error  propagation  are 


v'  _  A 

6  T.  =  >  - J-  6  N.  +  >  - 

J  “■**  mm  1 


1=3  *Ni 


0  T. 
1=J  l 


6  T. 

1 


N 


be. 


N 


5c.  ,y_j_  6n.  +  y  — - 

J  XTM  1  X  *■ 


be.  __ 

J—  6  T . 


1=3  *Ni 


bT. 
i=3  i 


j  =  K  N  (18) 


(Note  that  these  equations  are  triangular  because  the  value  of  a  transverse 
variable  at  z  does  not  depend  on  the  form  of  any  radial  function  for  r  <  z.  ) 
The  appropriate  partial  derivative  matrix  for  error  propagation  is  then 
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The  inverse  equations,  Eq.  (17),  are  not  known  explicitly  and  thus  the  partial 
derivatives  of  Eq.  (19)  cannot  be  calculated.  The  partial  de rivatives  bN/bT, 
bN/bc,  bT/bT  and  b^/bc  can,  however,  be  computed  with  the  radiation 
calculation  routines  of  EMABIC.  The  derivatives  are  computed  numerically. 
For  example,  the  array  of  derivatives  bN,/bT.  (i  =  1,  N)  are,  in  principle, 
computed  from 

bN.  n*  -  n: 

_ l  _  i _ l  (20 


where 


N+  =  N.  [T.,  T. , . , . . .  ,  (1+s)  T., 
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The  profiles  and  (i  =  1,  N)  are  computed  with  the  Tj  and 
obtained  from  the  inversion,  e  is  a  small  parameter  (e  =  0.01  is  generally 
used).  The  process  is  repeated  for  all  j  to  get  the  full  matrix  of  BN^/ BT^. 

The  three  remaining  derivative  arrays  are  computed  in  the  same  manner. 

In  practice,  the  derivatives  are  computed  in  a  slightly  different  manner. 

The  derivates  considered  so  far  refer  to  changes  of  transverse  variables  at 
transverse  zone  centers  with  respect  to  variations  at  radial  zone  centers,  or, 
vice  versa.  The  computation  routines  of  EMABIC  do  not  operate  on  variables 
defined  at  zone  centers,  but  rather  at  zone  boundaries.  Thus,  Eqs.  (20)  and  (21) 
are  used  to  compute  the  four  partial  derivatives  corresponding  to  the  corners 
of  the  intersection  of  the  ith  transverse  and  jth  radial  zone,  and  the  derivative, 
Eq.  (20)  approximated  as  the  average  of  these  four  values.  The  reason  for  oper¬ 
ating  on  zone  centers  instead  of  zone  boundaries  is  discussed  below. 

From  the  derivatives  just  discussed,  the  matrix 
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is  constructed.  Its  inverse  is  the  desired  propagation  matrix  A,  that  is, 

A  a  B'1  (23) 

The  inversion  of  B  is  carried  out  with  a  standard  Gaussian  elimination 
algorithm  using  maximal  pivot  strategy.^ 

The  objection  to  formulating  the  error  propagation  formalism  on  zone 
boundaries  is  that  the  resulting  derivative  matrix  corresponding  to  (22)  is 
singular  and  cannot  be  inverted.  The  singularity  results  because  partial 
derivatives  evaluated  at  i  =  j=  N+l  are  identically  zero. 

The  variance -covariance  matrices  for  propagation  through  the  inversion 
matrix  A  were  treated  in  Section  2.  2,  but  separately.  In  order  to  be  consistent 
with  A,  which  accounts  for  the  coupling  of  N  and  t  ,  the  appropriate  error 
matrix  is  constructed  from  Mj-j  an(j  M—  as  the  2N  x  2N  matrix 


The  assignment  of  zero  to  the  off-diagonal  quadrants  of  this  matrix  reflects  the 
usual  experimental  condition  that  the  radiance  and  absorptance  data  are  measured 
independently  from  each  other.  Thus  there  is  no  correlation  between  N.  and  T. 

f  •  •  1  J 

even  for  l  =  j. 

6.  B.  Carnahan,  H.  A.  Luther  and  J.  O.  Wilkes,  Applied  Numerical 
Methods,  John  Wiley  and  Sons  Inc.,  N.  Y.  ,  1969,  p.  282. 


Finally,  the  variance-covariance  matrix  for  the  recovered  temperature 
and  concentration  profiles  is 


(25) 


The  first  N  elements  on  the  diagonal  of  Mrpc  are  the  recovered  temperature 
variances,  and  the  last  N  elements  are  the  concentration  variances.  Explicitly, 
these  are 
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3.  EXAMPLE  APPLICATION 


An  empirical  analysis  was  made  in  a  previous  study ^  of  error 
propagation  for  temperature,  H^O  concentration  and  CO^  concentration 
retrieval  from  E/A  data  obtained  on  a  Transtage  engine.  The  present  error 
propagation  formalism  is  applied  here  to  temperature  and  CO.,  concentration 
retrieval  for  the  same  conditions  as  the  previous  analysis.  From  a  number  of 
successful  inversions  of  the  Transtage  data,  the  temperature  and  CC>2  concen¬ 
tration  profiles  of  Fig.  1  were  obtained.  In  the  previous  empirical  analysis, 
these  profiles,  along  with  the  pressure  profile  of  Fig.  2  were  used  to  generate  the 
transverse  E/A  profiles  of  Fig.  3.  Random  errors  were  then  superimposed  on 

the  E/A  data  at  the  M=201  equally  spaced  points  between  z  =  0  and  z  =  R=67cm. 

2 

The  rms  magnitude  of  the  errors  was  3%  of  peak  value  (cr=0.  0048  W/cm  -sr-pm) 
for  radiance  and  5%  of  peak  value  (o=0. 0093)  for  absorptance.  These  profiles 
are  shown  in  Fig.  3  and  essentially  model  the  original  unsmoothed  Transtage  data. 
The  noisy  E/A  profiles  were  then  smoothed  and  inverted  to  retrieve  the  radial 
Tc  profiles.  The  process  of  error  superposition,  smoothing  and  inversion  was 
carried  out  five  times.  For  all  five  cases,  smoothing  parameter  values  of 
vrr  =  675  and  v—  =  2700  were  used.  The  resultant  five  sets  of  retrieved  T  c 

N  T 

profiles  were  then  used  to  compute  the  standard  deviations  of  T  and  c  as  a 
function  of  r.  The  results  are  shown  in  Figs.  4  and  5. 

In  the  present  analysis,  no  smoothing  or  inversion  is  necessary.  We  can 
proceed  directly  from  the  radial  profiles  of  Figs.  1  and  2  (needed  to  compute  the 
partial  derivatives  bN/bT,  bN/bc,  br/bTand  b^/bc)  and  the  data  for  M,  R,  o(N.), 
ct(Tj),  Yj^  and  y—  already  given.  The  resultant  standard  deviations  cr(T.)  and 
cr(c.)  i=l,  N  for  N=5,  10  and  20  are  shown  in  Figs.  4  and  5. 

The  results  are  essentially  in  agreement  with  those  obtained  in  the  previous 
analysis.  The  present  analysis  yields  a  values  that  are  ~2  times  larger  than 
the  empirical  results.  This  difference  may  be  due  in  part  either  to  the  limited 
sample  (5)  used  to  compute  the  empirical  results  or  to  the  fact  that  the  empirical 
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Fig. 


Transtage  Temperature  and  Gas  Concentration  ] 
(a)  Temperature;  (b)  Concentrations. 


2 


1 


Fig.  5 


procedure  allows  for  nonlinear  error  propagation  while  the  present  analysis 
employs  a  linear  approximation. 

The  variation  of  the  results  with  N  shows  only  an  increasing  structure 
with  increasing  N  but  no  consistent  trend  to  change  the  overall  magnitude  of 
the  estimated  errors. 
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APPENDIX 


ERROR  PROPAGATION  CONTROL  CARDS  FOR  EMABIC 

The  incorporation  of  error  propagation  into  EMABIC  requires  a  slight 
modification  and  enlargement  of  the  instructions  for  data  preparation  given  in 
Chapter  IV  of  Ref.  1. 

The  first  change  allows  the  input  of  experimental  error  data  along  with  the 
transverse  data  itself.  In  Fig.  6  (p.  65)  of  Ref.  1,  (i)  and  o— (i)  should  be 

entered  in  columns  31-40  and  41-50,  respectively,  on  the  same  cards  that  enter 
N(i)  and  T(i).  The  legend  of  Fig.  6  should  include  the  following: 

a=r:(i)  =  estimated  rms  error  of  N(i)  TT  ..  .  .  .  ... 

N'  '  '  Units  consistent  with  N(i) 

_  and  ~  (i) 

a— (i)  =  estimated  rms  error  of  t  (i) 

The  second  change  concerns  the  List  Control  Card.  In  the  section  on  card 
type  10  (p.  64),  the  last  sentence  should  read,  "If  the  variable  JKLIST  (format 
A10)  on  the  LIST  card  has  the  value  PRINT,  the  emission  and  absorption  functions 
J(r)  and  K  (r)  are  also  listed  for  each  iteration.  " 

The  additions  are  two  new  control  cards  named  ERROR  and  NOERROR. 

The  following  addition  should  be  made  to  Fig.  3  (p.  58) 


and  the  following  section  added  after  the  discussion  of  the  Execution  Control 
Card  (pp.  66-67). 
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13.  Error  Control  Cards  Two  cards  .vith  the  names  ERROR  and 
NOERROR  control  error  propagation  calculations.  An  ERROR  card 
indicates  that  an  error  propagation  analysis  should  be  made  in  either  the 
Abel  or  iterative  Abel  inversion  mode.  If  TYPE  (format  TJC)  has  the  integer 
value  0  or  blank,  the  error  analysis  will  be  performed  immediately  following 
a  successful  inversion  run  in  which  E/A  data  have  been  inputed.  An  error 
analysis  does  not  require  an  inversion.  If  TYPE  has  the  value  1,  an  error 
analysis  will  be  performed  using  the  error  data  read  in  from  the  transverse 
data  deck  and  the  pTc  data  read  in  from  a  profile  generation  mode  deck.  A 
SPECIE  card  must  be  used  in  this  mode. 

The  TYPE  =  1  mode  allows  an  error  analysis  to  be  made  for  hypothetical 
experiments.  For  example,  if  first-order  estimates  of  radial  profiles  are 
already  known,  and  likely  experimental  errors  are  known  from  previous 
experience  (or  are  estimated),  it  can  be  determined  whether  or  not  an 
experimental  program  is  likely  to  improve  on  the  first-order  estimates  of 
the  radial  profiles. 

The  NOERROR  card  turns  off  any  error  analysis  for  subsequent  runs. 
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